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I. INTRODUCTION 



Many attempts 0-3] have been made to account for the unexpected behavior of the recent expansion of the universe. 
In a previous paper [ 1 1, we showed that non-perturbative quantum effects in the vacuum may account for the SNe-Ia 
data suggesting a recent acceleration of the expansion of the universe. We considered a free quantized scalar field 
(with possible coupling to the scalar curvature, R). The propagator of the field included an infinite sum of all terms 
having at least one factor of R |t],|8| . The effective action was given in Q , which can be referred to by the reader for a 
more complete exposition of the theory. We found an approximate solution to the effective Einstein gravitational field 
equations. In that solution, there was a sudden transition (at redshift z of order 1) from an earlier matter dominated 
stage of the expansion to a mildly inflating de Sitter expansion. This sudden transition model required the universe 
to be spatially open in order to fit the cosmological data. 

Here we find an improved analytic solution to the effective Einstein equations and determine its consequences. 
Surprisingly, we find that the analytic solution permits the SNe-Ia data to be fit with a spatially-flat model having 
a reasonable matter density and age. The spatial flatness makes a significant difference || in attempting to fit the 
spectrum of small angular scale fluctuations of the cosmic microwave background radiation (CMBR). It appears that 
the existing data favors a spatially-flat cosmological model. Spatial flatness is also favored in models of the universe 
having an early inflationary stage. 

In our present model, there is only one free parameter, namely, the ratio of the mass parameter rn to the present 
value of the Hubble constant, Hq. It is convenient to express this parameter in the form 

m h = m/h, (1) 

where 



h = iV(100km/(sMpc)), (2) 

with Hq expressed in km/(s Mpc). We find that a good fit to the SNe-Ia data is obtained when m,h is in the range 
6.40 x 10 -33 eV < rrih < 7.25 x 10 - 33 eV. For this parameter range, the ratio of the total matter density (including 
dark matter) to critical density, flo, is found to be in the range 0.58 > Hq > 0.15, and the age to of the universe is 
found to be in the range 8.IO/1 -1 Gyr < to < Y1.2YT 1 Gyr. Also, we find that our model gives reasonable abundances 
for light elements formed during big-bang nucleosynthesis. In our model, no special coincidence at the present time 
is necessary to explain why the energy density of matter is of the same order as the vacuum energy density. 

The organization of this paper is as follows. In Section II, we give the effective Einstein equation of the model, and 
summarize the solution of Ref. |Q. In Section III, we derive an analytic solution for the scale factor of the model. In 
Section IV, we compare the model to SNe-Ia data and obtain ranges for mh, ^0 and ^o- I n Section V, we discuss the 
implications of the model for big-bang nucleosynthesis and for the spectrum of CMBR fluctuations. In Section VI, 
we show it is probable (independent of fitting the observations) that, at the present time, the vacuum energy density 
in our model is comparable to the matter density. Finally, our conclusions are given in Section VII. 



II. OUR MODEL 

In Ref. Q, we consider a free, massive quantized scalar field of inverse Compton wavelength (or mass) m, and 
curvature coupling £. The effective action for gravity coupled to such a field is obtained by integrating out the 
vacuum fluctuations of the field . This effective action is the simplest one that gives the standard trace anomaly 
in the massless-conformally-coupled limit, and contains the explicit sum (in arbitrary dimensions) of all terms in the 
propagator having at least one factor of the scalar curvature, R. 

The trace of the Einstein equations obtained by variation of this effective action with respect to the metric tensor 
take the following form in a Friedmann- Robertson- Walker (FRW) spacetime (in units such that c = 1): 

R+I±-4A = J^-{(rn 2 +lR)hi \ l+lRm' 2 ' 



m2 ~ (R M^ + ^ir-aosor.^H, (3) 



m 2 + £R V 2 m 

where T c / is the trace of the stress tensor of a classical, perfect fluid component containing mixed matter and radiation, 
A„ is the cosmological constant, k q = (167rG) _1 (G is Newton's constant), £ = £—1/6, and v is a quantity that vanishes 
in de Sitter space: 
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(4) 



Here, the metric is 

rf s 2 =-d( 2 +fl (f) 2 ^ +r wj, (5) 

where a(i) is the scale factor and k = 0, 1, —1 give spatially flat, closed and open universes respectively. 

The right-hand-side (RHS) of Eq. (||) is proportional to the trace of the quantum contribution to the stress tensor. 
In Eq. (g) , we have assumed that terms in the quantum stress tensor that depend on derivatives of the curvature are 
negligible. This assumption is consistent with the solution that we obtain. 

In 141, we found that Eq. (§) admits a solution in which a matter-dominated FRW universe transits to a stage of 
the expansion in which the scalar curvature is nearly constant. This led us to construct an approximate cosmological 
model which makes a sudden transition to a de Sitter universe out of a matter-dominated one. In that approximation, 
we needed a spatially-open cosmology in order to fit high-redshift SNe-Ia data, and to get a reasonable value of the 
present matter density. 

In this paper, we eliminate the assumption of a sudden transition to a de Sitter expansion, by finding an analytic 
solution to the effective Einstein equations of our model. Consistency with cosmological observations requires that this 
improved analytic solution has zero (or nearly zero) spatial curvature. This solution consists of a matter-dominated 
universe smoothly, with continous first and second derivatives of the scale factor (i.e., C ), joined to a constant scalar 
curvature expansion that asymptotes to a de Sitter universe at late times. The spatial flatness is consistent with 
estimates of the first acoustic peak of small angular scale fluctuations in the CMBR, which suggest that the universe 
is spatially flat (although this issue awaits an ultimate decision at the time of writing, see Ref. @ for details). 

We now turn to the derivation of a smooth solution for the cosmological scale factor, a(t). 



III. ANALYTIC SOLUTION 

In the following analysis, as in [Q, we take A G = 0. For £ < 0, and with a low value of m (~ 10 -33 eV), we find in 
H that the evolution of the FRW universe is essentially unaffected by the quantum contributions at early times. A 
perturbative analysis in h shows that quantum contributions to the stress tensor begin to have a significant effect at 
a time t j , when the classical scalar curvature has decreased to a value given roughly by 

Rcl{tj ) = *£il=m 2 , (6) 

where 

m 2 = m 2 /(-£). (7) 

This effect is shown to occur well into the matter-dominated stage of the evolution (therefore, for t > tj, T c i ~ —p m ). 
We further argue that, after the scalar curvature reaches the value R c i (i.e. for t > tj), it stays essentially constant 
during the later evolution of the universe. More precisely, for all times t > tj, we find that the scalar curvature has 
the form 

R(t) = m 2 (l-e(t)), (8) 

where 

<t) = ^(-S+(S 2 +4f3r 1 ) 1/2 ), (9) 



and 



pm(t) --i _ r 

I nn - 

2160£ 



Here. 
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mpi being the Planck mass. As shown in it is straightforward to verify that e(t) is of order r for t > tj and of 
order y/r at t = tj. For the low value of r under consideration, it follows that e(t) <C 1 for all t>tj. Thus the scalar 
curvature stays nearly constant at the value m 2 for t > tj. 

Despite the constancy of the scalar curvature, the quantum contribution to the stress tensor increases dramatically 
from the time tj to the present time to. To see this, consider the trace of the quantum stress-tensor, T q , defined by 
T q /(2n ) = -(RHS of (ff). It follows from (g) that (with A Q = 0), 

T q = -2k q R - T c ( 

= -2 Ko M 2 + p m + 0(e), (12) 
where the second approximate equality holds for t >tj. Assuming (as we find) that the transition time tj occurs at 



a redshift Zj of order 1, and given that p m (tj) ~ 2n m 2 (from Eq. (§D), we obtain the matter density at the present 
time, p m (t ), as 

Pm(to) = -T cl {t ) = (a(^)/a(t )) 3 p m {tj) ~ (1 + z,) _3 2K ra 2 

~ K m 2 /4 + 0(e). (13) 

Thus, according to Eq. (|l2|), T q grows from T q (tj) ~ —2K m 2 + 2K m 2 ~ at time tj to T q (t a ) ~ — 2K Q m 2 + 
(1/4)k to 2 ~ — (7/4)k to 2 ~ 7T c i(t ) at the present time to- Thus, we find that the quantum contribution to the 
stress tensor grows from a negligible value at a redshift of order 1 to a value exceeding the classical contribution at 
the present time. 

To obtain the scale factor a(t), we consider a spatially flat (k — 0) cosmology, and find the constant scalar 
curvature solution for t > tj that joins to the usual matter-dominated solution at time tj, with continous first and 

matter-dominated universe, is 



second derivatives. The time tj obtained from Eq. (o), using R cl (t) = (4/3)t~ 2 for t < tj in the spatially flat 



tj = (2/v / 3)m _1 . (14) 
For such a universe the Hubble constant H(t) = (2/3)t~ 1 , has the value at tj given by 

H{tj) =m/V3. (15) 
For t > tj, we now obtain the unique solution to the constant scalar curvature equation 

R = <i(- + K) =™\ (16) 
\a a z ) 

that satisfies 

H(t 1 ) = d(t ] )/a(t J ) = m/V3. (17) 

Continuity of the scalar curvature at tj ensures that the second derivative of the scale factor is also continuous. This 
solution is 



a(t) = a(tj) \ I sinh ( —f= — a J / sinh [ — — a ) , t > tj, (18) 



V3 J I V3 



where 



a= H -tanlT 1 ~ 0.117. (19) 

We have verified that this solution also satisfies the remaining Einstein equations up to terms of order e. According 
to Eq. (|l8|), the expansion approaches a de Sitter expansion at late times (i.e., tfn S> 1). Furthermore, it has the 
property that the deceleration parameter go = —cia/a 2 is negative (i.e., the universe is accelerating) for all t > t a , 
where 
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t a = VSmr 1 (a + tanh -1 ^ -1 / 2 )) ~ 1.73m- 1 ~ 1.50 f,-. (20) 
Also, the solution joins in a smooth (C 2 ) manner to the usual spatially flat matter-dominated solution for t < tj, i.e. 

a(t) = a(tj) (t/tj) 2/3 = (\/3mi/ 2 ) 2/3 , t < tj. (21) 

This solution, given by Eqs. (|l8|), ( [i"9| ) and (pl|), depends on only one parameters to. It is more accurate than 
the approximation of a sudden transition to a dc Sitter expansion, that we used in |4|. Our analytic solution takes 
a long time (of order m -1 ) to approach a de Sitter expansion. This is the reason that the spatially flat model gives 
agreement with observations (as we show below). In the approximation of a sudden transition to de Sitter, consistency 
with observations required a negative spatial curvature. We further note that our analytic solution is not the same 
as that of a mixed matter-cosmological constant model, for which a(t) oc (sinh(\/3A</2)) 2 ^ 3 - 

Further improvements in the accuracy of our solution can be made by an iterative procedure in which the present 
solution is used to evaluate e(t) in Eqs. (||) and (|To|), and then Eq. (|l^) is solved with m 2 replaced by the RHS of 
Eq. (§). 

IV. COMPARISON WITH OBSERVATIONS 

For comparison with observations of high-redshift Type la supernovae (SNe-Ia), we calculate the luminosity 
distance-redshift relation for the model defined by Eqs. (|l^) and (^l]), and, from it, the difference in apparent 
and absolute magnitudes as a function of the redshift z of a source. This difference is given by 

m- M = 51og 10 dx +25, (22) 

where g?l is the luminosity distance to the source in Mpc, defined as |ll]] 

d L = (l + z)a r 1 , (23) 

where ao is the present scale factor, and n is the comoving coordinate distance from a source at redshift z to a 
detector at redshift 0. For a spatially flat FRW universe, ri is given by 

i r dz' , s 

"—7. mp-y m 

where H (z) is the Hubble constant as a function of z, and the speed of light c is now shown explicitly. It is convenient 
to introduce the parameter h = i?o/(100 km/(s Mpc)), where Hq is the present value of the Hubble constant in km/ (s 
Mpc). It is also useful to work with the rescaled mass parameter mh of Eq. (Q). We note that to/j, being an inverse 
length, can be measured in units of Mpc -1 . 

For our model, we define zj as the redshift at time tj. Then we obtain, using Eqs. (pi) and (p4|), 



hri<(z) = a 1 Vl2m h 1 x J 
where ri < (z) denotes r\(z) for z < zj, and 



dx 



(25) 



X — sinh [ ° - — a I . (26) 



V x/3 

For z > Zj, we use Eqs. ( pl| ) and (24) to obtain 

hr x> {z) = h ri< (z 3 ) + 2%/3a^ 1 TO- 1 (l + Z] f' 2 ((1 + z^) -1 / 2 - (1 + z) -1 / 2 ) . (27) 

Eq. ( |23| ) gives the luminosity distance for this model, dLi(z), defined by 

d L i(z) = (1 + z)a ri < (z), z < Zj 

= (l + z)a r 1> (z), z>z r (28) 
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After substituting for ri<(z) and r% > (z) in the above equations, hdLi(z) does not depend on ao (for a spatially flat 
universe). The three parameters that occur in hd^i are hto, Zj and writ- However, we may differentiate Eq. ( |l8| ) to 
obtain 

cm , / ctm 



H(f) = 7!I C0 'H7f- t 'J' (29) 

The above equation, evaluated at the present cosmic time to, can be rewritten in the form 

ht = (3.26 x 10 6 yr/Mpc)™,; 1 (tanh" 1 ((865.4 Mpc)m h ) + a) , (30) 
which gives ht in years. Also, Eq. yields 



1 + zj — a(to)/a(tj) — J sinh ( - — °__ - — a ] /sinh ( — — a ) . (31) 



V3 // \3 

Eq. ( p0| ) shows that hto is a function of alone. Therefore, the redshift Zj as given by Eq. ( |3l| ) is also a function of 
to?, alone. This implies that the function hdt,x(z) depends on a single parameter, nih- The present ratio of the matter 
density to critical density, Slo, is also a function of the same parameter. To see this, we use p m oc a~ 3 . Then one has 

8ttG , , 8ttG , . / HUAV fa(t,j)\ 3 , . 



Continuity with the spatially flat matter-dominated universe at t = tj requires that 
{%TtGp m {t 3 ))/{3H{tj) 2 ) = 1. Therefore 

O = (H^/Ho) 2 (afe)/a ) 3 

= (2.996 x 10 6 Mpc 2 )m^ ^sinh ^ c/l ^ n,t ^ a ) / sinh a ) ) ' ( 33 ) 



where we have used Eqs. (15) and (|18j) in arriving at the last equality. Note that the quantity hto appearing in the 
above equation is itself a function of (see Eq. (|30|)). Therefore flo is also a function of the parameter m^, which 
is the only adjustable parameter in our model. 

If the time tj is less than the present age to, the monotonic behavior of Eq. ( p9| ) implies that H(oo) < Ho < H(tj). 
With H(oo) and H(tj) obtained from Eq.([29|), the previous inequality gives 

to > V3H /c= 5.78 x l(T 4 /i Mpc" 1 , (34) 
m < Vl2iV c = 116 x 10~ 3 h Mpc" 1 . (35) 

We therefore find that the rescaled mass parameter mh is constrained by the model to lie in the range 

5.78 x lO^Mpc^ 1 <m h < 1.16 x lO^Mpc^ 1 . (36) 

The above range on mh can be also expressed in electron-volts (eV), as 

3.69 x 1(T 33 cV < m h < 7.39 x 1CT 33 eV. (37) 

We emphasize that the above range of values of to/j is a consequence of our model, if to > tj , independent of any fit 
to cosmological observations. 

Plots of flo and hto vs. rrih, using Eqs. (^) and (|30| ) respectively, are given in Fig. |l| for the range of to/j values 
of Eq. (37) above. For the same range of to/j values, Fig. || is a plot of flo vs. hto- 

We now use the luminosity distance of Eq. (^8[) to fit the SNe-Ia data. To do so, it is convenient to normalize the 
difference to — M of Eq. ( |22| ) to its value in an open universe with f^o = 0.2. We define 

where 
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d L 2(n , z) = 2if - 1 cfV 2 («o* + (n - 2) (Vl + O * - l) ) 

= h- 1 (5995.8 MpcjQ- 2 (ct z + (fi - 2) (V 1 + - l)) • (39) 

Since hdi,i(z) is a function of m/j and z, and ft,di2(0.2, z) is a function of z alone, it follows that A(m — M) is a 
function of rrih and z. 

Fig. |^ is a plot of A(m — M) vs. z, along with a plot of SNe-Ia data acquired from Ref. 0. A good fit is obtained 
by any curve that lies in between the two dashed curves shown, where curve (a) has rrih = 6.40 x 10~ 33 eV, and curve 
(b) has rrih = 7.25 x 10~ 33 eV. This gives a value of £Iq in the range 

0.15 < Q < 0.58, (40) 

and a value of hto in the range 

8.10 Gyr < ht < 12.2 Gyr. (41) 

Estimates of h give 0.55 <h< 0.75 §. For h = 0.65, Eq. © gives a range for the age of the universe, namely 

12.5 Gyr < t < 18.8 Gyr. (42) 

A representative solid curve (c) has the value rrih = 6.93 x 10~ 33 eV, which gives fi = 0.346 and t — 14.8 Gyr 
(with h — 0.65). A plot of the scale factor a(t) for this curve is shown in Fig. [|. 



V. RECOMBINATION, NUCLEOSYNTHESIS AND CMBR FLUCTUATIONS 

In this section, we first calculate the time of recombination in our model, compare it to the recombination time in 
a standard open universe, and discuss the implications for big-bang nucleosynthesis. We then show that the apparent 
angular size of CMBR fluctuations in our model is somewhat smaller than in a spatially flat fio = 1 model, although 
consistent with available data. 

In our model, we find that recombination occurs in the matter-dominated era. The redshift at recombination is 
given by the following expression, discussed in Ref. |12|, 

z r = 1048 (l + 0.0124(fVi 2 )- ' 738 ) (l + gi (n h 2 ) 92 ) , (43) 

where 

9l = 0.0783(n b h 2 r - 238 (1 + 39.5(fi^ 2 ) - 763 ) _1 , (44) 

g 2 = 0.560 (1 + 21.1(f2 b /i 2 ) 1 - 81 )~ 1 , (45) 

and SI;, is the ratio of baryon density and critical density at the present time. For the ranges 0.0025 < f^/i 2 < 0.25, 
0.15 <tt < 0.60 and 0.55 <h< 0.75 §, we find, 

1055.73 <z r < 1301.80. (46) 

The redshift at matter-radiation equality, z eq , is given by 

i + z eq = n h 2 /(n r h 2 ), (47) 

where tt r is the ratio of radiation energy density to critical density at the present time. With a CMBR temperature 
of 2.726 K at the present time |]l3| , we obtain fl r h 2 = 4.16 x 10~ 5 . Using the same range of values as above for f2o 
and h, we find 

1090.14 < z eq < 8115.90. (48) 

Therefore, in our model, recombination occurs in the matter-dominated era. 

We now invert Eqs. (JTI^) and ( ^l| ) to obtain a range for the time of recombination, t r , corresponding to the range 
of z r in Eq. ([i"6|). A lower bound on t r is obtained with the higher value of z r and the lowest value of consistent 
with a fit to the SNe-Ia data, i.e., m,h = 6.40 x 10 -33 eV, and an upper bound is obtained with the lower value of z r 
and the highest value of rrih consistent with the same data, i.e., rrih — 7.25 x 10 - 33 eV. We find 
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1.82 x 10 5 years < ht r < 4.97 x 10 5 years. (49) 

We compare the above range with the corresponding range of ht r in a standard, spatially open matter-dominated 
model, with the same ranges of fto, h and ili,h 2 (and therefore the same range of z r ). The scale factor for an open 
matter-dominated model pT{ | leads to the parametrized equations 

2(1 - Op) , , . n _i 
1 + z r = — (coshi/v — 1) , 

"0 

ht r = (4.89 x 10 9 yr) _"° o)3/2 (sinh W - A) , (50) 

tftr being a dimensionless parameter. The above equations give the following range for ht r in an open matter-dominated 
universe, 

1.79 x 10 5 years < ht r < 4.89 x 10 5 years. (51) 

Comparison of Eqs. (|i^) and ( |5l|) reveals that the range of recombination times for the two models differ by only 
about 1.7 percent. Furthermore, the scale factors prior to recombination are almost identical in both models (the 
effect of spatial curvature in the open model is negligible for times prior to the recombination time). Therefore, it 
is expected that nuclear reaction calculations in big-bang nucleosynthesis would give nearly identical results for the 
abundances of light elements in both models. Since these calculations are known to give results in agreement with 
observation for the standard open model, nucleosynthesis calculations in our model would also give results that agree 
with observations, for the same values of fio, h and flbh 2 . 

We now turn to the calculation of the apparent angular size of a fluctuation of a fixed proper size D at the surface 
of last scattering. This apparent angular size in our model will be compared to the apparent angular size of the same 
fluctuation in a spatially flat FRW model having Qq = 1, since the latter model is known to yield a CMBR fluctuation 
spectrum consistent with data. 

For simplicity, we take the representative values = 6.93 x 10 -33 eV, h = 0.65, and Vl^h 2 — 0.025 in the discussion 
that follows. The chosen value of corresponds to the intermediate curve (c) of Fig. | and gives Qq = 0.346. Eq. 
( [43| ) then gives the redshift at last scattering, as z r — 1089. 

The apparent angular size 9 of a fluctuation of fixed proper size D at last scattering is given by 

9 = D/(a(t r )n) = D(l + z r ) 2 /d L1 (z r ), (52) 
where dLi(z r ) is found from Eq. (p8[), which also gives 

Mil (1089) = 1-017 x 10 7 Mpc. (53) 

Therefore, with h = 0.65, 

9 = (7.594 x 10~ 2 Mpc^D (54) 

in our model. 

In a spatially flat f^o = 1 model, taking h — 0.65 and Q^li 2 = 0.025, Eq. ( [43| ) gives z r — 1105. 
The luminosity distance for a spatially flat Q = 1 cosmology is given by 

hd L {z) = (5995.8 Mpc) (1 + z) (l - (1 + z)~ 1/2 ) , (55) 

which gives 

hd L (1105) = 6.433 x 10 6 Mpc. (56) 

Therefore, with h = 0.65, and using Eq. (|52|), we obtain 

9 = (1.236 x 10" 1 Mpc" 1 )!?. (57) 

Comparison of Eqs. ( |54"|) and (|57j ) reveals that the apparent angular size of a fluctuation of a given proper size at 
last scattering is about 1.63 times less in our model than in a spatially flat H,q — 1 model. If the first acoustic peak in 
the small angular scale CMBR spectrum arises from a fluctuation of fixed proper size in all models, then the above 
result implies that, in our model, this peak would be shifted to a higher mode number relative to that in a fio = 1 
spatially flat model. In the latter model, the first peak is known to occur at mode number I ~ 200, and therefore 
it would occur at I ~ 326 in our model. These two possibilities are both consistent with the existing data on small 
angular scale CMBR fluctuations (see, for example, B and references therein). 
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VI. THE QUESTION OF FINE-TUNING 



Recent observations of Type-la supernovae evidently indicate that the universe has been in an accelerating phase 
from a redshift of order 1 up to the present time, which implies that the contribution of matter to the total energy 
density is of the same order of magnitude as the contribution of vacuum energy density at the present time (in 
spatially flat models, this means that f2o/(l — fio) is of order 1). Why should this be the case? In mixed matter 
and cosmological constant models, this question requires an explanation of why the cosmological constant must be 
fine-tuned to a precise, non-zero value. As pointed out in Q, our model is relatively insensitive to the value of the 
cosmological constant term. However, in our model, it would appear that, within its allowed range of Eq. (|37|), the 
parameter rrih must be finely tuned to give values for Zj and — £Iq) that are within an order of magnitude of 1. 

It should be noted that the allowed range of values of rrih given by Eq. ([57j) does not by itself constrain Zj and f2 . 
The lowest allowed value of rrih gives Zj = and Q,q = 1, and the highest allowed value gives Zj = oo and £lo = 0. 

Here, we argue that, in our model, values within an oder of magnitude of 1 for Zj and ^o/(l — ^o) are more likely 
than other values. The argument that follows rests on two assumptions: (i) to > tj> and (ii) all values of rrih within 
the allowed range given by Eq. ( p7| ) have equal a priori probability. Assumption (i) implies the range of rrih values 
of Eq. (|37]), and assumption (ii) is reasonable in lieu of a detailed fundamental theory that predicts the value of rrih- 

Given these two assumptions, one may compute the probability distributions for Zj and f^o in a straightforward 
manner, since both quantities are functions of rrih alone. By assumption (ii), the probability distribution function for 
to/i, P(m^), has the form 

P(m h ) = P , 3.69 x l(T 33 eV < m h < 7.39 x l(T 33 eV, 

= 0, otherwise. (58) 

Normalization of P(rrih) yields P = 2.70 x 10 32 / eV. We then obtain the probability distribution functions P(zj) and 
P(£!o/(l — f^o)) for Zj and ^o/(l — ^o) respectively, as 

p{z ' )=p ° n drt: 1 (59) 

P(O /(l - n )) = Pa/ I d(no/ . (1 = " 0)) |, (60) 

dm h 

where Zj(rrih) and fio( m /i) are given by Eqs. ( pl| ) and ( |33"| ) respectively. We now compute the probability that Zj and 
f2o/(l — ^o) lie between 0.1 and 10 (i.e., they are within an order of magnitude of 1). We find 

P[0.1 < z 3 < 10] = / dzj P(zj) = 0.851, 
Jo.i 

p[o.i < n /(i - n„) < io] = J d (v^A p (t%^) = °- 619 - ( 61 ) 

It is therefore more likely for Zj and — Oq) to be within an order of magnitude of 1 rather than outside that 

range, assuming that all allowed values of rrih have equal a priori probability. 



VII. CONCLUSIONS 



We have shown that vacuum effects of a free scalar field of very low mass can account for the observed acceleration 
(i.e., the SNe-Ia data), while at the same time predicting reasonable values for the age of the universe and the total 
matter density. Evidently, our model also predicts reasonable light element abundances, and as a consequence of 
spatial flatness is in agreement with current data on small angular scale CMBR fluctuations. 

Better SNc-Ia data would be able to distinguish between our model and mixed matter-cosmological constant models, 
as well as quintessence models ||. The curves of A(m — M ) vs z are different in our model from those of the other 
models. Future observations of small angular size CMBR fluctuations may also distinguish between these models. 

We emphasize that our model is based on a free renormalizable quantum field and does not require that we live at a 
very special time in the evolution of the universe. We also note that a graviton field of very low mass may give rise to 
vacuum effects similar to those of the scalar field we considered here. In contrast, a scalar particle of very high mass 
with similar characteristics to the present one, may contribute to a stage of early inflation, with reheating and exit 
from inflation caused by particle production from the vacuum and possibly from other inflationary exit mechanisms. 
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FIG. 1. Plots of O and hto versus m h , with the range of m h values 3.69 x 10 33 eV < m h < 7.40 x 10 33 eV. 
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FIG. 2. A plot of Q versus ht , with the range of m h values 3.69 x 10 33 eV < m h < 7.40 x 10 33 eV. 
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FIG. 3. A plot of the difference between apparent and absolute magnitudes, as functions of redshift z, normalized to an open 
universe with Qq = 0.2 and zero cosmological constant. The points with vertical error bars represent SNe-Ia data obtained from 
Ref. [5]. The two dashed curves represent the values (a) mt = 6.40 x 10 -33 eV (lower dashed curve), and (b) ran = 7.25 x 10 -33 
eV (upper dashed curve). The solid curve represents the intermediate value (c) m/, = 6.93 x 10~ eV. 
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